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Abstract 

The Lee- Yang theorem for the zeroes of the partition function is not strictly apphcable to quan- 
tum systems because the zeroes are defined in units of the fugacity e^'^'^ , and the Euchdean-time 
lattice spacing At can be divergent in the infrared (IR). We recently presented analytic arguments 
describing how a new space-Euclidean time zeroes expansion can be defined, which reproduces 
Lee and Yang's scaling but avoids the unresolved branch points associated with the breaking of 
nonlocal symmetries such as parity. We now present a first numerical analysis for this new zeros 
approach for a quantum spin chain system. We use our scheme to quantify the renormalization 
group flow of the physical lattice couplings to the IR fixed point of this system. We argue that the 
generic Finite-Size Scaling (ESS) function of our scheme is identically the entanglement entropy of 
the lattice partition function and, therefore, that we are able to directly extract the central charge, 
c, of the quantum spin chain system using conformal predictions for the scaling of the entanglement 
entropy. 
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I. INTRODUCTION 



The quantum spin chain system has been important for the development of much of the 
modern understanding of the role and nature of quantum fluctuations in low-dimensional 
quantum spin systems. The three important results for this system were derived in the 1980's 
and 1990's from treatments parameterised in terms of local fluctuations of the vacuum angle 
6*, at finite values of 6', where 9 is defined as the prefactor of the topological term, 

Q = J dtjdx n.{dxn x dtu) (1) 
with n a state vector for the spins in the chain, and t and x the imaginary-time and spatial 
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coordinates of the model [ij [2] • Firstly, it was found that the infrared (IR) fixed point of 
the quantum spin chain system at 6' = vr is directly related to the unstable renormalization 
group flow of the SU(2) Wess-Zumino-Witten model, such that the central charge of the 
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quantum spin chain system at^ = 7risc=l [JJ. Secondly, it was identified in [3| that the 
quantum spin chain system is asymptotically free and, thirdly, it was established in J] that 
the action of the quantum spin chain can be described effectively by the Sine-Gordon model 
in the vicinity of the IR fixed point at = tt. These three results indicate that conformal 
symmetries are important for classifying the quantum spin chain system, and also, that a 
nonperturbative description of the system would be needed in order to treat the quantum 
fluctuations of quantum spin chain systems at arbitrary values of 6. 

Formally, the procedure of renormalization defines how the scale of the fluctuations of 
a given system are related to the physical couplings. However, no procedure has yet been 
defined to renormalize the fluctuations in 6 of the quantum spin chain at arbitrary 6. Al- 



though it was established in ^\ that perturbative renormalization is a well defined program 
for the effective action given by the Sine-Gordon model at the IR fixed point, it is not clear 
the extent to which this renormalization program can be generalised to arbitrary values of 6. 
In fact, it has been argued more recently in [5| that the nonlocal fluctuations of 9 are equally 



as important to the quantum spin chain system as the loca^ 



a 



fluctuations considered previ- 



we have recently argued in 



ously. Since the quantum spin chain is asymptotically free 
6| that the full renormalization program for the quantum spin chain at arbitrary 6 requires 
a nonperturbative renormalization scheme. The purpose of this article is to establish such a 
nonperturbative renormalization prescription for the quantum spin chain and, consequently, 
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to quantify the link between the physical couplings of a quantum spin chain experimental 
system and the underlying conformal symmetries. 

One way to treat the nonperturbative dynamics of the quantum spin chain is via a nu- 
merical lattice simulation using the Quantum Monte Carlo (QMC) method. This approach 
defines a realisation of the nonperturbative vacuum of the quantum spin chain that is given 
by evaluating the lattice partition function of the system via importance sampling. The 
approach generates a finite statistical ensemble to represent the lattice partition function, 
which consists of a set of different configurations of the lattice generated at a fixed value of 
the inverse temperature, /3, the interaction coupling between the spins, J, and the length 
of the quantum spin chain, L. There are two immediate problems with trying to define the 
nonperturbative renormalization group flow of the quantum spin chain via the QMC nu- 
merical procedure. Firstly, the quantum spin chain systems generated by the QMC method 
are restricted to finite values of /3, whereas the IR fixed point of the system, identified in 
hJl, is at /? = oo. Secondly, it is not possible to give an area preserving definition for ^ as a 
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function of the lattice volume [2]P^PJJ] and, consequently, to define how d fluctuates as a 
function of lattice system size. However, we have recently resolved both of these issues alge- 
braically in and by defining a new Wick rotation procedure for the local transfer 
matrix elements of the QMC method, in order to define an exact mapping between (3 and 9. 
Although there is no critical point in the quantum spin chain system at finite temperature, 
we are able to define the line of constant physics in the (3 — J plane that flows into the 
quantum critical point from this scheme. 

Our new Wick rotation procedure is essentially ageneralisation of the properties of par- 
tition function zeros, as we discussed in detail in [7|. The reason for this is that both of 
these schemes are based on a mathematical procedure called meromorphic continuation. 
Partition function zeros are conventionally defined through an exact polynomial expansion 
of the lattice partition function defined over the local degrees of freedom of a lattice system 
12] , and the zeroes of the finite polynomial expansion defined for a finite volume lattice sys- 
tem uniquely define the singularities of the lattice partition function of that finite system. 
However, there is also a nontrivial relationship between the zeroes of a finite lattice system 
and the zeroes of that same system treated in the infinite lattice volume limit [l^. This 
relationship is defined formally by a meromorphic mapping. 

To define our new Wick rotation procedure, we have used the formal definition of a mero- 
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morphic mapping to treat a more generalised polynomial expansion than the conventional 
local polynomial expansion, which includes the nonlocal degrees of freedom of the lattice 
model of a quantum spin system in addition to the local degrees of freedom. As we have rig- 
orously proven in we can, therefore, now use this new polynomial expansion to quantify 
the singularities of the partition function of a lattice system associated with the nonlocal 
fluctuations of 9, as well as the local fluctuations of 9. This determination defines our new 
nonperturbative renormalization program for fluctuations in 6*, since it completely defines 
the fixed points of the nonperturbative vacuum in 9 in the mapping to the infinite lattice 
volume limit. The numerical analysis we present is, therefore, simply a repetition of the 
conventional Finite-Size Scaling (FSS) analyses given in 27|- 3^ for the conventional parti- 
tion function zeros of classical spin systems, but now applied to our new form of nonlocal 
polynomial expansion. The only sig nificant difference that arises between this new zeros 
analysis and the analyses in j27|-j32| is that the new generic FSS function of the zeros has 
the form of the entanglement entropy of the system, as we discuss in Section 7. This means 
that we are then able to compare the results of our new zeros analysis with conformal results 
on the FSS of the entanglement entropy of the quantum spin chain, and from this we deduce 
the central charge of our lattice ensembles. 

One practical point remains, which is that although our Wick rotation procedure defines 
the meromorphic continuation properties of 9 on the lattice, we would still like to be able to 
vary the fixed value of 9 at which our lattice ensembles are generated in order to generate 
lattice ensembles at different values of 9. The reason for this is that this will allow us to 
assess the sensitivity of our new renormalization approach to perturbations in 9. Although 
it would be possible for us to generate lattice ensembles directly at the point 6' = tt by 
simulating models consisting entirely of half-integer spins with the same spin magnitudes, 
it would not be possible to vary 9 by varying any of the lattice coupling parameters in these 
models. We, therefore, choose to simulate the model which consists of an antiferromagnetic 
(AFM) periodic mixed-spin chain of the form, 1 — 1 — 3/2 — 3/2, where 1 and 3/2 represent 



the spin magnitudes, and the period of the alternation is four lattice sites 13|] We can 



vary the fixed value of 9 at which a lattice ensemble is generated in this model by varying 
the ratio of the nearest-neighbour spin interaction couplings within the periodic cell. The 
Hamiltonian of this model is given by. 
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TT _ \^ J , T o(l) q(3/2) . <j(3/2) <j(3/2) . <j(3/2) ,^(1) 

-H - 2^ '^l,l'^4j •'^4i+l + -'l,3/2'54j+l-'^4i+2 + '^3/2,3/2'54j+2-'^4j+3 + l-^J 

j=0 

where j is the lattice site index of the cell along the spin chain, Ja^b is the nearest neighbour 
spin interaction coupling between the neighbouring spins of magnitude a and b, and S''-'^^ is 
the spin operator of magnitude a. For our numerical simulations of the model we keep the 
ratio between like-like and like-dislike spin vectors fixed such that Ji^i = ^3/2,3/2 = J and 
Ji.3/2 = J3/2.1 = oiJ ■ We organise the article as follows. In Section 2 we review the definition 
of the transfer matrix of the continuous-time QMC method, and in Section 3 define how 
this can be used to give an exact polynomial expansion for the lattice partition function 
of the quantum spin chain. In Section 4 we discuss the symmetry properties of our new 
polynomial expansion, and how the zeros of the polynomial relate to the zero temperature 
fixed points of the quantum spin chain through Wick rotation. In Section 5 we identify the 
generic FSS behaviour of our new zeros density function, which we identify as the entropy of 
the QMC method. In Section 6 we present the numerical FSS analysis of an AFM periodic 
mixed-spin model of the form 1 — 1 — 3/2 — 3/2. Finally, in Section 8, we identify that the 
entropy that we have measured is identically the entanglement entropy of the system, which 
follows from the bipartite definition of the partition function obtained from Trotter-Suzuki 
decomposition, and we compare our numerical FSS result with conformal predictions for 
the scaling of the entanglement entropy at the 6 = n fixed point of the quantum spin chain 
system. 



II. THE TRANSFER MATRIX 



A. The Loop-Cluster Lattice Partition Function 

There are several closely related loop-cluster QMC methods for the numerical study of 
low- dimensional lattice quantum spin systems, and a comprehensive review of all of these 



schemes is given in 16[. Trotter-Suzuki decomposition is used to define the generic partition 
function for a number of these QMC methods, and in particular, the method we have used 
to generate our numerical data in this article. In this decomposition the spin operators, 
S*, of a Z) dimensional lattice quantum spin system are placed on a il* + 1 dimensional 
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space- Euclidean time lattice as classical spins. The spin Hamiltonian, Ti = J^^SiSi^i, for 
a system defined with an AFM nearest-neighbour interaction is then split into a sublattice 
definition on the space-Euclidean time lattice, where the spin site index i of the Hamiltonian, 
Hi = JSiSi^i, is defined only on the even-indexed spin sites, with J is the AFM nearest- 
neighbour spin interaction coupling. Including different spin vectors does little to change the 
basic form of this decomposition, since higher order spins can be implemented locally through 
the inclusion of a permutation factor 9j. The lattice partition function of the generic D 
dimensional quantum spin system with an AFM nearest-neighbour spin interaction is defined 
via Trotter-Suzuki decomposition through the following trace. 



JJexp(-Ari/, + C(Ar2 



(3) 



i=l 



where (3 is the inverse temperature. At is the Euclidean-time lattice spacing, i is the 
lattice site index, and L is the number of sites on the lattice 16|. This recovers the quantum 



spin operators of the D-dimensional model from the classical spins of the D + 1-dimensional 
model. For our analysis of the quantum spin chain D = 1. This choice of sublattice 
Hamiltonian used for the Trotter-Suzuki decomposition of a quantum spin system with an 
AFM nearest-neighbour spin interaction means that the lattice partition function can be 
expressed via a basis of two-by-two site blocks known as plaquettes [l^ , 



Z = Tre^^= H (5,,i5,+i,i|e-^^«|5,,i+i5, 



(4) 



where (.) denotes the transfer matrix for one of the elementary plaquettes defined in i^t, 
where i is the spatial lattice site index (in the D direction), t is the index of Euclidean-time 
sites on the lattice (in the -|-1 direction), and T is the total number of sites in the Euclidean 
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time (+1) direction [16|. 

The transfer matrix (.) of an 0(3) spin system with a nearest neighbour exchange is of 
the form, 



/l 







1(1 + e-^^^ 



1(1 



^JAt 







-J At 







e-^^n 



(5) 



1/ 
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In this definition the diagonal elements describe the possible Euclidean-time evolution 
of a spin defined on the spatial sites, and the off-diagonal elements describe what must 
subsequently happen to the spins on the nearest-neighbour spatial sites in order to locally 
conserve total spin via an antiferromagnetic exchange. The form of the transfer matrix that 
is given in is essentially the same as that used for the Density Matrix Renormalization 
Group, but the important difference for the QMC method is the way in which the above 
Trotter-Suzuki formalism can be used to define loop-clusters for the Monte Carlo procedure. 
The general approach of the QMC method is closely related to classical loop-cluster methods, 



as is discussed in detail in |l6|. 



B. Continuous-time 

The general aim of applying the Monte Carlo scheme to lattice models is to generate a 
sequence of lattice configurations via a dynamical Markov process. An importance sampling 
decision separates each of the lattice configurations in this sequence (or Markov chain), and 
the set of all of these lattice configurations then defines a statistically meaningful ensemble 
and the lattice partition. The important consequence of this generic Monte Carlo approach 
for the evaluation of (jlj) is that there are then in fact three discrete indices for the lattice 
partition function evaluated via the generic Monte Carlo method. The first of these indices 
is the spatial index, i, the second is the Euclidean-time index, t, and the third (which is 
usually dropped from notation) is the configuration index, M, which indicates the relative 
position of a lattice configuration within the Markov chain. 

New lattice configurations are generated for the importance sampling step of the Monte 
Carlo process by making small random changes to the state of the spins on the lattice, and 
in the classical loop-cluster scheme these changes are made by changing the state of the spins 
on the lattice within a closed loop. The point of using this loop-cluster definition for the 
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importance sampling of classical spin systems, as was treated rigorously in [171 [18|, is that 
the sum over discrete steps in the probabilistic Markov chain, M, can then essentially be 
interchanged with the sum over discrete loop-clusters in the definition of the lattice partition 
function, which leads to more efficient numerical sampling. As was noticed in and 16|. 
however, when this loop-cluster definition of the new lattice configurations is applied to the 
Trotter-Suzuki formalism in (jlj), this definition has the advantage that closed loops on the 
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D + 1 dimensional lattice in (jl]) automatically preserve the definition of the trace. Therefore, 
as is used explicitly in [l^, using this approach it is possible to interchange, t, with, M, 
and to define the lattice partition via the with respect to a trace operation defined over the 
length of the Markov chain, rather than with respect to Euclidean-time. A formal definition 
of the continuous-time properties of the general continuous-time Monte Carlo method is 
given in 19|, and it is explained in detail how, t, and, M, can be interchanged in a general 



Monte Carlo method by using continuous-time. 

The QMC method we employ for our zeroes procedure is known as the continuous-time 
QMC method, and it corresponds to an algorithm where, t, is interchanged with, M, in 
the definition of (jl]). We define the lattice partition function of the continuous-time QMC 



method from 



151 ] by the following product of local transfer matrices. 



Z = Tie 



I3H 



n 



(6) 



o\ 

A, 1 - A, 
1 - A, A, 
\0 1 / 

where i is the spatial lattice site index, and Aj are the local transition probabilities 
corresponding to the entries of the transfer matrix. To be clear, these transition probabilities, 
Aj, are taken from the final figure in 15|, but there is no explicit expression in 15| that defines 
the lattice partition function as a product of these matrix elements. Therefore, to obtain 
(ini) we have substituted these transition probabilities from the final figure into the transfer 
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matrix definition in (jlj) (which is also from [15||). Our reason for using a more general form of 
the transfer matrix than ([5]) is that we want to distinguish the transfer matrix elements that 
are realised numerically for each lattice site, z, since these values are explicitly evaluated at 
each step of the Markov process and are used to build up a realistic probability distribution 
for the spins. 

The most obvious difference from the way the above lattice partition function is written 
from (jlj), however, is that the trace over the spatial and Euclidean-time directions of the 
D + 1 dimensional lattice in (jlj) have been replaced in (j6j) by a trace over just the spatial 
index, i. In fact, to write (j6j), the trace over t = {1,T} has been replaced in (j6j) by a 
trace over just one site, M = {1,2}, through the exchange t ^ M. To write the lattice 
partition function in this form we have used the continuous-time properties of the general 
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continuous-time Monte Carlo method defined in [19|. We have noticed, firstly, that both 
the number of discrete sites in Euclidean-time, T, and the number of discrete steps in the 
Markov chain, M, are arbitrarily defined for the partition function in (jlj). Secondly, for a 
Markov process, the importance sampling process represents a probabilistic decision that is 
continuously defined in the interval (0, 1]. Therefore, because of the more rigorous functional 
definition of the probabilistic interval that the steps in the Markov chain represent, there is 
no restriction on the number of steps in the Markov chain, and it is possible to represent 
the Markov chain by only one step. This is also the rigorous basis of dropping the usual 
statistical configuration index on lattice partition functions evaluated via the Monte Carlo 
method, [l9]. To be clear, this is not the usual way that the boundary conditions of the 



continuous-time QMC method in [15| are presented in review articles 20] 21]. However, the 
usual tacit assumption that is presented: that M commutes with t, is only strictly true in 
the large volume limit. In a rigorous sense, the continuous-time method is a nonequilibrium 
Monte Carlo method 13], as we have presented it here, where the noncommutativity of M 
and t can be realised in the numerical simulation of a finite system. 

The value of expression ([6]) for evaluating our new scheme is twofold. Firstly, it is straight- 
forward to define an exact polynomial expansion for the lattice partition function in ([6]) 
(defined over i), as we will now show. Secondly, the I At I are defined to be a continuous 
functions over (0, 1] in the continuous-time method 15|] 19|], which allows us to define our 
exact Wick rotation procedure for the lattice system. 



III. LATTICE PARTITION FUNCTION POLYNOMIAL EXPANSION 

To form a polynomial expansion over the spatial site index i for the lattice partition 
function in ([6]), we first introduce, P, which is a 4z x 4z block diagonal matrix, where the 
blocks are the local transfer matrices, (.), that define the matrices in the product in ([6]) and 
all other entries in P are zero. 



(.)m/ 

{.)k+l,M 









{.)l,m ) 



(7) 



9 



An exact polynomial expansion for the lattice partition function, in /3, is given by finding 
the characteristic equation of the eigenvalue problem det(P — /3) =0. 



There is a precedent for using numerical 



the zeros of the lattice partition function 



attice simulation data in this way to calculate 



22| 23|. There are difficulties for practical imple- 



menting these schemes because the evaluation of the characteristic polynomial coefficients, 
Cj, from numerical input from standard numerical recursion can be subject to roundoff errors. 
This a general problem for all exact diagonalization and Density Matrix Renormalization 
Group analysis that involve diagonalization, and often this roundoff forms the dominant 



systematic error, as is discussed at length in 2^ 2^ 26 1. This source of roundoff error is 
both a function of the number of sites in a lattice model and also the condition number of 
the system that is numerically realised. Often, this roundoff restricts the range of lattice 
volumes that can be treated via an approach involving diagonalization. To address this 



problem in our analysis, we have implemented the tested techniques in 22| designed to min- 
imise roundoff. These techniques involve using Newton's relations to define the operation 
d = Cj[Aj], in combination with a series of logarithmic modifications to the standard nu- 
merical operations. We refer the interested reader to these studies 221] 23|] (and references 
therein) for further details of the numerical stability of these techniques, but comment that 
the size of the coefficients that we have evaluated for the quantum spin chain system mean 
that the roundoff problems from the routines themselves are negligible as a function of lattice 
volume since they have been used to treat successfully treat polynomials of (9(2500). 

The application of Newton's relations is a standard linear algebra procedure for calculat- 
ing characteristic polynomial coefficients. The procedure consists of evaluating Cj = Cj[Aj] 
recursively from powers of the trace of P, as defined in 22|. We introduce our Wick rotation 
at this point in the analysis, mapping (3 ^ i(3 through the following expression for the trace 
of the local transfer matrices in (El), 



Tr (.) = det 



1 \i 
\i 1 



(8) 



This Wick rotation definition, following 6|], [71] and [8[, changes the boundary conditions 
of the lattice partition function from (0, (3] —>■ [— vr, tt] . The trace of P is then easily calculated 
from the product of these local transfer matrix traces via. 
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TrP = ^(1 - A^) (9) 



i=l 

The characteristic polynomial equation can then evaluated via Newton's method to find 
the polynomial coefficients in the Wick rotated basis, and the lattice partition function 
polynomial expansion can be written as, 

L+l n-1 

Z = J2{c^)f3^' : nc„ + TrP" + ^CfcTrP""^' = (10) 

i=l k=l 

where i the spatial lattice site index, (3 is the inverse temperature, q are the polynomial 
expansion coefficients, n = 1, L, and /c is a dummy index. 

To summarise the steps we have taken to define this exact polynomial expansion for the 
lattice partition function of the quantum spin chain: we have first taken the definition of 
the continuous-time lattice partition function in [l^ and formulated this expression (|6]) as 
eigenvalue problem in (3. We have then Wick rotated the local transfer matrix elements 
following the algebraic procedure we have defined in j^, [3] and ^, by defining the trace 
operation in ([8]), to map [3 ^ i[3 and re-define the boundary values of the lattice partition 
function from Euclidean to imaginary time. Finally, we have then proposed evaluating the 
characteristic equation of this eigenvalue problem using Newton's relations and the numerical 



recursion modifications defined in 



221 and 



231 . This will then give a semi-analytic exact 



polynomial expansion, defined over the lattice spatial site index which can be evaluated 
from numerical input in the form of the local transition probabilities Aj realised in the 
importance sampling of the continuous-time simulation of the quantum spin chain system. 



IV. PROPERTIES OF THE POLYNOMIAL ZEROES 



In the previous section, we have described our analytic procedure for converting the 
numerical transfer matrix entries of the continuous-time QMC method into an exact semi- 
analytic polynomial expansion for the lattice partition function of a quantum spin system. 



We claim that the zeros of the lattice partition function in 
properties to the partition function zeroes defined in [271 ]- 
algebraic Wick rotation procedure we have defined in 



0]^ should share similar scaling 
321 ] because we have used the 
and Q] to define the polynomial. 
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One of the most important properties of the zeroes defined in [27[-[32| comes from an 



argument presented in the Appendix of 



271]. This states that if the local elements used 



to define the characteristic polynomial coefficients are bounded in magnitude, then the 
magnitude of the coefficients is necessarily bounded, as is the magnitude of the zeroes. This 
result is important to the FSS analysis of partition function zeroes because it means that 
the zeroes are then constrained to lie on a locus in the complex expansion parameter plane, 
which simplifies the definition of the FSS function The generic FSS function of the 

zeroes is given formally in the next section. However, it can be noted that for our new 
expansion in ( fTOl) a similar bound exists because our polynomial coefficients are defined 
from the local transition probabilities Aj, which are bounded in the interval (0, 1]. 

For our Wick rotated procedure, the meaning of this locus constraint has a special con- 
nection with the way that 0(3) spin operators of the quantum spin chain are represented in 
the continuous-time method. Changing the continuous-time representation of 0(3) spins via 
the Wick rotation procedure in (IHl) means that the S'^i-component of the 0(3) spin vectors is 
represented by the spin vectors defined on each spatial lattice site i, and the Euclidean-time 
interval, (0,/9], represents the projection of the spins onto the plane of the Sx-Sy compo- 
nents. Therefore, the lattice nearest-neighbour spin interaction coupling J only couples to 
the S';;-component of spin of the new Wick rotated operators, and the coupling between the 

and Sy components is implicitly normalised to unity. The bounds that define the locus of 
our new expansion are, therefore, set by the physical surface traced out by the spin vectors. 
In an isotropic system, where the interaction couplings to Sx, Sy and Sz are all unity, we 
can deduce that the zeros of (ITOl) should lie on a sphere because the spin vectors have a 
47r-rotational symmetry, ie. the locus is a sphere. However, any deviations of J from unity 
will break the isotropy of this surface traced out by the spin vectors. We can, therefore, 
include the analytic factor J into the exact polynomial in (fTOl) in order to normalise the 
units of the surface for our analysis, 

L+l 

z = Y,{cHPJr- (11) 

i=l 

This then allows us to relate the zeroes locus to the lattice volume (defined in physical 
units), and to make the connection between the critical scaling of the lattice partition 
function in 9 and the physical units of the lattice ensemble. 
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V. THE ZEROES DENSITY 



The most important feature we have claimed for our expansion in (fTOj) is that it shares 



the same FSS properties as the partition function zeroes analyses defined in [27|]- 32|. We, 
therefore, first review the FSS function arguments for these classical spin systems, and then 
present a new FSS definition for our new expansion that we show is equivalent. 



A. The Zeroes Density of the Ising Model 

The FSS function of the generic partition function zeroes analysis is defined by the density 
of the partition zeros along the locus that they form in the c omp lex plane of the expansion 



parameter 



271]- 32|. For the Ising model system considered in 27|, the zeroes density, A(0, t) 



is defined through a saddlepoint equation for the free energy, F, of the following form 



-f3h 



X{(f),t) logle" -e"^]d(j) 



(12) 



where, h, is the applied external field in the Ising model, t is the reduced temperature, is 
the angle that a zero on the locus subtends to the real e'^-axis, and the expansion parameter 
for the partition function zeroes treatment is the fugacity variable e^. The zeros density, 
X{(j),t), is also subject to the boundary conditions. 



A(0,t) = A(-0,t) ; 



1 



A(0,t) d0 = - 



(13) 



The only stable solution for the free energy in f|T2|) occurs at = and this indicates 
that the critical points of the Ising model considered in 12| must lie on the real e'^-axis. 
This critical point is the saddlepoint solution for the free energy, found in the asymptotic 
limit where the zeroes density becomes a continuous function. What is expected on a finite 
lattice is that none of the zeroes will lie on the real e''-axis, but that the zeroes density is 
still an asymptotically convergent function towards the critical point. This assertion is used 
to argue that the zeroes density function is a valid FSS function in 12|. Note, however, that 
there is no rigorous statement in this definition that proves that (fT2l) is the free energy of 
a lattice partition function. The statement to justify the partition function zeros method 



given m 



12| is simply that the thermodynamic properties of a lattice system can be treated 



exactly if the free energy of a lattice partition function is written in the form of ( |T2l) . 
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B. The Zeroes Density of the Quantum Spin model 



Unfortunately, we cannot make a simple substitution to map our new expansion variables 
into f[T21) in order to define the zeroes density of our new formulation. The reason for this 
is that the boundary conditions of our new formulation are branch points, because of the 
way we have defined the Wick rotation which maps [0,(3] [— 7r,7r]. We start by noticing 
that the Aj in (6) have an explicit dependence on the couplings jS and J. This dependence 
arises through the bare parameter dependencies of the plaquette definition in ([6]), which are 
used to construct the Aj through importance sampling [l^[2^. Each of these numerically- 
realised local transfer matrix entries represents a transition probability, therefore, for these 
transition probabilities to be well-defined it follows that the integral over all values of the 
bare parameters must be one. This motivates a normalizing condition, analogous to f|T3l) . 



oo 



A(/3,J)>0; / \{p,J)dJ = l. (14) 







In our notation the Aj are the local transfer matrix elements of a finite, discrete lattice 
system, whereas the function X{(3, J) is the transition probability corresponding to these 
matrix elements defined for a continuous system, as a function of the continuous variable 
J. Note that the /3-dependence of the local transfer matrix elements Aj is already made 
continuous by the definition of the continuous-time method in [lol {20]. 

A rigorous treatment of the existence of this continuous-J limit of the transition proba- 



bilities is given by applying the method of ^-function renormalization 34], which has been 
developed for this lattice model in [3]. We define an analogous saddlepoint equation to the 
one in ( |T2l) for the polynomial expansion in (ITOl) . by defining the analytic properties of the 
support of this lattice model. We do this via quantifying the logarithm of the determinant 
of the operator P in ([7]), where. 



detP = JjAi. (15) 



i=l 



and Aj are the eigenvalues of P, Aj = 2Aj — 1. To quantify this logarithm (whilst resolving 
any undetermined branch points in the continuous-system limit, L 00) via the (^-function 
renormalization prescription, it is necessary to analytic continue the following function, 
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(16) 



1=1 



from large s to s = 0. Introducing the ^-function, Cp{s) = '^i=i ^7^^ logarithm of 
the determinant of P is completely defined by the first derivative of this C- function at s = 0. 
Writing this ^-function in the usual form of a Mellin transform 3J| 35|| 36|| , 



Cp{s) = Tr (e-^^) J-^ d/ (17) 



r(.) 

fOO 



X{f3, J) dJ (18) 







/■CXD 

IndetP = - / A(/3, J) ln(J) dJ (19) 
Jo 

Physically, this expression for the logarithm of the determinant of P in f|T9l) is related to 
the entropy of the quantum spin chain, and is the analogue of (fT2|) for this quantum system. 
The logarithm of the determinant of P is equivalently the logarithm of the partition function 
Z, from (6), and since there is an implicit integration over (3 defined for the continuous-time 
method integrating over P yields a prefactor of /S, 

/ In detP d(3 = (3\n detP = H (20) 
Jo 

where H is of the form of an entropy. However, unlike the classical spin system in (fT2l) . 
the total entropy of the quantum spin system is not fully defined from this expression in (ITQll . 
and H in (20) only represents the local limit of the entropy of the quantum spin system. In 
physical terms, this quantum system differs from the classical system in that it is defined via 
a time-ordered path integral over /3, which in analytical terms, is expressed by the fact that 
not only is the first moment of f|T9l) finite, but also a number of the higher moments of the 
C-function are finite. We made an explicit discussion of this problem in 33||, commenting 
that the free energy maxima and saddlepoint solutions of a quantum spin system need not 
be coincident in the general case. 

To determine the maxima of the logarithm of the determinant of P, and hence the 
maxima of the support of the system, we have to evaluate the (^-function prescription at a 
different point on the fundamental strip, at s = 1. This is essentially to include the nonlocal 



contributions to the path integral within a finite analytic domain of J 

15 



35| . where. 



\s=i= A(/3, J) W{J) d/, J = 0...L (21) 
The full expression to be analytically continued to determine this support, from (16), is, 

POO 

H= A(/3, J) lnA(/3, J) dJ (22) 
Jo 

which corresponds to the more familiar Von Neumann form of the entropy of a quantum 



spin system [37|. This expression is not as analytically well defined as the local expression 
for the entropy in (fTOjl . because the branch points in the path integral over /3 are unresolved 
at s = 1. Hence, in the general case, the Von Neumann entropy can be a divergent function 



However, since we know from the ^-function prescription in (1191) that it is possible to 
write down a local limit for the entropy, it is sufficient (to determine the properties of the 
support that we are interested in) to show that the maximum value of this local limit is 
unique. In physical terms, this means that we can treat the important case of the quantum 
system having a degenerate groundstate, by identifying the maxima of the support of the 
system but without specifying which of these degenerate vacuua we are considering. 

In order to establish that both our new definition of X{(3, J) is a valid FSS function we 
must, therefore, establish the location of the maxima of the Von Neumann entropy in (22) 
is asymptotically convergent towards a particular point in J as a function of discrete J. To 



do this we follow the analytic treatment given in [35|. The total entropy of the quantum 
spin system is constructed from the method of Lagrange undetermined multipliers from the 
(local) moments in (21). Firstly, X{(3, J) is approximated by a series that contains a finite 
number of degrees of freedom M, 

A(AJ)~exp|^-^7,(/5)J^j (23) 
Then, the moments of (22) are defined in order to explore the branch structure of J, 



/ij- = / In^J A(/5, J) dJ (24) 
Jo 

The maxima of (22) is then asymptotically convergent to one fixed value of J, as a 



function of M, provided that the following functional is convex. 
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M 



r(7i(/5), ... , 7Af (/?)) = Yl ^^(f^) /^J- + /^o In 

3=0 



- rexp(-^7,(/?)ln^^/) dJ 



What this means is that the Von Neumann entropy in (22) then has a well defined 
asymptotic limit as a function of any finite local degree of freedom M, such as lattice 
system size, since the normalizing conditions in (14) guarantee the convexity of the above 



functional 



35|. 



Both H and A(/3, J) are therefore valid FSS functions for our numerical lattice calculation, 
and on a practical level, the continuous probability density function A(/9, J) can be expressed 
by the number density of the local transition probabilities Aj. We can therefore identify (19) 
and (14) as our quantum system analogues of the defining relation for the zeros density of 
the classical spin system in fll3p . and its normalizing conditions in f[T21) . Our definition for 
the quantum analogue of the zeros density is the continuous probability density function 
X{f3, J). From (19) the zeroes of our new expansion follow exactly the same form of FSS as 



the zeros of classical systems treated in 



271 ]- 32], and so we can choose a simple difference 



relation to define A(/3, J) for our finite lattice systems 311] 



A(0) = A 



1 



(25) 



2 ; N{(f)k+i - <pk) 

where k is the sequential index assigned to the zeros along the locus they form in the 
complex- (/? J) ^ plane, and is the angle subtended from the real-(/5J)^ axis to a given zero. 
Since A(/3, J) is a well defined FSS scaling function, the scaling of the zeros density describes 
the renormalization group scale transformations of the lattice system. 



A(/5, J, L) = L^XipL, JL) 
where c is the scale transformation parameter. As is usual for FSS functions 



(26) 



39||40| 



m 



the vicinity of a second order critical point (|26|) . this leads to us approximating that. 



lim A(J) = ril-aJ... 

J,I3~*0 



(27) 



where a is an arbitrary constant. As the /3-dependence of A(/?, J) is continuous, by 
construction, we can also make the substitution e**^ = J to define. 
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lnA(z0) = cln(J) — aJ ... ; lnA(z0) ~ c/ 
Similarly, the first order scaling contribution is of the form, 



(28) 



^^ffi ^ I A(i0) = const. (29) 

This form of this function follows from the residue theorem applied to the point where 
the zeroes density develops its maxima, which is a nonanalytic point in the infinite lattice 
volume limit. We label this point, Jq, and some general contour around this point in the 
complex (/9J)^-plane, C = e*'^ + Jq. Therefore, 

— — = Am^TiR, (30) 



ic J ~ Jo Jo 
VI. NUMERICAL RESULTS AND ANALYSIS 

We have generated lattice ensembles the AFM periodic mixed spin chain of the form, 
1 — 1 — 3/2 — 3/2, where 1 and 3/2 represent the spin magnitudes and the period of the 
alternation is four lattice sites, by using the continuous-time method |l^ 1^. We consider 
this model defined with two independent nearest neighbour spin couplings, defined via (1). 
We keep the ratio between the like-like and like-dislike neighbouring spin pairs fixed in our 
analysis and consider the analytic continuation of the lattice partition function in j3 and 
J. Separate treatments of this model via the Nonlinear Sigma Model (NLSM) and Density 
Matrix Renormalization Group (DMRG) methods have indicated that the groundstate of 
this system should have three distinct topological regions, with a gap f ul g roundstate, sepa- 



411 ] 42| corresponding to 



41| 42|. However, neither of these analyses 



rated by two gapless quantum critical points at zero temperature 
the vacuum angles 6 = n and 6 = dir. In this model 6 can be varied by varying the ratio 
of the nearest neighbour spin interaction couplings within the periodic cell. If one of these 
couplings is fixed to unity the two quantum phase transitions are found at values of the 
second coupling of = 0.483 and Jc = 1.314, 
has been able to quantify the renormalization group flow of the physical couplings to the 
6 = 7T IR fixed point. 

For our numerical analysis we have generated a range of lattice data sets each consisting 
of 100, 000 decorrelated lattice configurations at fixed /?, J and chain length L using the 
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continuous-time method. We have then evaluated the polynomial expansion coefficients 
Ci defined in f fTOl) by using the procedure defined in Section 3, using the local transition 
probabilities realised in the importance sampling of the continuous-time simulation of the 
quantum spin chain system as our numerical input. We have used a standard rootfinder to 
then calculate the zeros of these polynomials for each of the lattice ensembles. These zeroes 
are plotted for two of the ensembles in Fig.l. The general picture is consistent with our 
expectation: namely, that the zeroes are constrained to lie on a regular locus in the complex 
expansion parameter plane, as we argued in Section 4. The zeroes in Fig.l are plotted for 
lattice ensembles generated at two different values of J. The first of these values, Jc = 0.483, 
corresponds to the region around the quantum phase transition at ^ = vr whereas the second 
value, Jc = 1.314, corresponds to the region around the quantum phase transition at ^ = 37r. 

A. Critical values of /3J 

We find in Fig.l that both of the two lattice ensembles not only contain information 
about the phase transition points that the lattice ensembles are generated in the vicinity of, 
but also the second transition that the ensembles are a considerable distance from in J in 
the phase space. There is a critical point (indicated by the point at which the zeroes density 
(12^ sharply increases) in both ensembles at both = and at = 0.9, where (f) is the angle 
subtended from the positive real-axis anticlockwise in the figure. One important feature of 
our semi-analytic polynomial expansion is that it enables us to analytic continue the lattice 
partition function as a function of (3 J. Since the zeroes density approaches a maxima at the 
real axis at, = 0, the radius of the zeroes locus, {pj)iocus, is related to the coordinates 
of the point on the line of constant physics in the (3 — J plane which passes through the 
quantum critical point at zero temperature. Naively, we can substitute the fixed value of 
one of the lattice couplings used to generate the ensemble (either [3 or J) into this value 
of the radius, {(3J)iocusi in order to define the coordinates of this point in the (3 — J plane. 
However, practically the lattice partition function has a finite cutoff in J beyond which the 
analytic continuation in J is undefined. We tabulate the radius of the loci in Table 1 for the 
different lattice ensembles we generated. We find that the radius of the locus is unity when 
the ensembles are generated at a J value less than J^. = 0.483, but that the radius of the loci 
is rescaled by a specific factor when the ensembles are generated at a J value greater than 
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FIG. 1: Partition function zeros of the quantum spin chain in the complex (/3J)^ plane, for two 
different lattice ensembles. The lattice size for both systems is L=64 and the inverse temperature 
j3 = 40, whilst one lattice ensemble is generated at J = 0.483 and the other at J = 1.314. 



Jc = 0.483. Furthermore, the product of this specific factor and the value of J at which 
the ensembles are generated is unity, as we tabulate in the final column of Table 1. This 
indicates that when we generate ensembles at a point in the (3 — J plane where the line of 
constant physics that passes through the quantum critical point is above the J cutoff, but, 
below the (3 cutoff, then the critical scaling of the partition function is not dimensionful in 
physical units. 

From the final column of Table 1 we find that the lines of constant physics which pass 
through the two quantum critical points are located at roughly = 0.483(1) and Jc = 
1.314(0) in the (3 — 3 plane. These values can be obtained independently from all the lattice 
ensembles we have generated, and are consistent with the critical values found in the DMRG 
and NLSM calculations in 



B. Critical Scaling Exponents 

We plot the logarithm of our new zeros density function (11) in Fig 2, which we have 
measured via the simple difference relation in (15). Two forms of scaling are evident, in 
line with our expectations. Firstly, in the vicinity of the two critical points at = and 
(f) = 0.9 the logarithm of the zeros density follows the second order scaling fit of fl28|) . We 
tabulate our fitting results to (128|1 in Table 2. Within the error bars of the analysis, the 
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L 


/? 


J 




J X {f3J)locus 


16 


100, 40 


0.483, 0.483 


1.080, 1.080 


0.582, 0.582 


32 


100, 40 


0.483, 0.483 


1.038, 1.038 


0.492, 0.492 


64 


100, 40 


0.483, 0.483 


1.018, 1.018 


0.487, 0.487 


128 


100, 40 


0.483, 0.483 


1.004, 1.004 


0.484, 0.484 


16 


100, 40 


1.314, 1.314 


0.702, 0.702 


1.101, 1.101 


32 


100, 40 


1.314, 1.314 


0.639, 0.638 


1.050, 1.050 


64 


100, 40 


1.314, 1.314 


0.608, 0.608 


1.025, 1.025 


128 


100, 40 


1.314, 1.314 


0.591, 0.591 


1.010, 1.010 



TABLE I: Dependence of the value of the zero nearest the real axis in the complex (/3J)^-plane, 
{f3J)iocusi on the numerical ensemble parameters: lattice size, L, inverse temperature P, and nearest 
neighbour coupling J. The couplings and results in the columns are matched pairwise across the 
table. The final column indicates the relation between the value of the couplings at the critical 
point, in physical units, and the lattice volume in physical units. 



fitted value of the exponent is consistent with a value of around one for all of the lattice 
ensembles. The stability of our determination of the exponent apparently improves with 
increasing spatial volume until we reach L = 128, where the errors then increase. The error 
bars on the measurements are in goo d ag reement with similar measurements made from the 
QMC method using spin-correlators 50|] SI] 52|. The error bars are also in good agreement 
with recent DMRG and exact diagonalization studies 45| 47| 49|]. Furthermore, the system 
size at which our analysis appears to be come unreliable directly corresponds to the range 
of volumes that can be treated in these recent studies, which indicates a potential common 
dominant source of systematic error in the numerical roundoff errors from diagonalization. 

It is interesting to note, however, that the lattice volume dependence of the numerical 
stability of our measurements appears to be correlated with the second form of scaling that 
we have measured in our data, which is of the form of an overall ^-independent first order 
contribution to the logarithm of the zeros density, (19). If we compare lattice ensembles 
generated at the same fixed /5 and J in Table 2, the scale of the Euclidean-time spacing 
fluctuation varies by about 1% when we go from a spatial volume of L = 64 to a volume of 
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L = 128 (in the final column). 

The generic origin of this first order contribution to the zeros density is slightly subtle. 
Naively, the physical volume of the lattice in the continuous-time method is fixed by choosing 
the fixed values of /5 and J at which the lattice ensemble is generated. Generating lattice 
data sets at different values of L, at fixed /3 and J, should therefore correspond to the same 
physical volume of system. However, because the role of Euclidean-time and the Markov 
chain are interchanged via the use of loop-cluster techniques in ([6]) it is the size of the loop- 
clusters that are fixed in physical units and not the Euclidean-time lattice spacing since 
the Monte Carlo measure is defined in the continuous-time approach by the free energy 
of the average loop-cluster and not by the free energy of the lattice system [15|[19|. Our 
measurements indicate that the expectation of the Euclidean-time lattice spacing changes 
disco ntinuously as we go from spatial volumes of L = 64 to volumes of L = 128 in our lattice 
ensembles. Although the physical lattice volume is fixed, in principle, because the lattice 
operators have been Wick-rotated the approach is sensitive to changes in the numerically 
realised IR cutoff of the ensemble in Euclidean-time. 

As we discussed in detail in jsl, what is expected for the renormalization group flow of 
the physical couplings of the quantum spin chains in the vicinity of the IR fixed point is 
that there is a Kosterlitz-Thouless transition line in the vicinity of the fixed point at zero 
temperature. Therefore, potentially what happens for the L = 128 volumes is that the line 
of constant physics that is determined from the analysis no longer passes directly through 
the zero temperature fixed point, for this ensemble, but instead it hits the zero temperature 
axis some distance from the fixed point. This would change the values of the fitted exponents 
for L = 128. However, whether the change in the numerical stability has a physical origin 
or not the additional first order contribution to the ensemble for L = 128 indicates that 
there is a change in the condition number of the system being modeled 2J] 25|] 26|] , and this 
directly effects the numerical stability of the analysis through roundoff. 



C. Numerical Stability 

Importantly, the FSS results obtained in this new analysis have been obtained from a 
spectral density definition of the zeroes density function, given in (16), where the scaling 
is defined as a function of the spectral index, k, rather than as a function of the lattice 
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FIG. 2: The density of the zeros in the complex- (/? J) ^ plane for an ensemble generated at 

L = 512, /3 = 40 and J = 0.483, plotted as a function of the angle, (/>, which measures the angle 
subtended by a zero on the zeroes locus from the real axis of the complex- (/J J) ^ plane (anticlockwise 
along the ellipse in Figure 1). Two sorts of scaling are evident in the plot: a constant behaviour 
corresponding to first order scaling p9l) where the zeroes density is flat, and a divergent second 
order scaling behaviour ()28p corresponding to the two second order transition points at (/> = and 
= 0.9 . 

volume, L. This means that we have obtained the critical exponents in Table 2 by using the 
information contained within single lattice ensembles generated at just one arbitrary value of 
the lattice size and couplings. We stress, we have not combined the information contained in 
lattice ensembles generated at several different volume sizes to obtain our exponent results, 
nor have we combined the information contained in lattice ensembles generated at several 
different couplings values, which would be usual for the FSS analysis of Monte Carlo data. 
We have in fact obtained 24 independent, consistent results for this scaling exponent from 
24 independently generated lattice ensembles of varying lattice sizes and with varying fixed 
coupling values. We have also performed a parallel QMC analysis of this model in jl^ . 
using conventional spin-correlator based measurement techniques, and have found that the 
computer time needed to perform this new measurement for one lattice ensemble at a fixed 
volume and fixed couplings is comparable to the computer time needed to generate one 
spin-correlator measurement for one lattice ensemble at a fixed volume and fixed couplings. 
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L 


/? 


J 


c(0 = 0) 


c((^ = 0.9) 




16 


40 


0.483 


1.008(0.078) 


1.008(0.078) 


0.629(0.217) 


16 


100 


0.483 


1.008(0.078) 


1.008(0.078) 


0.628(0.217) 


32 


40 


0.483 


1.002(0.027) 


1.002(0.027) 


0.627(0.149) 


32 


40 


1.314 


1.004(0.041) 


1.004(0.041) 


0.627(0.152) 


32 


100 


0.483 


1.002(0.027) 


1.002(0.027) 


0.627(0.149) 


32 


100 


1.314 


1.011(0.115) 


1.011(0.115) 


0.625(0.143) 


64 


40 


0.483 


1.002(0.017) 


1.002(0.017) 


0.626(0.107) 


64 


40 


1.314 


1.002(0.004) 


1.002(0.004) 


0.623(0.109) 


64 


100 


0.483 


1.002(0.016) 


1.002(0.016) 


0.626(0.106) 


64 


100 


1.314 


1.002(0.009) 


1.002(0.009) 


0.623(0.109) 


128 


40 


0.483 


1.216(0.218) 


1.216(0.218) 


0.616(0.088) 


128 


100 


0.483 


1.163(0.193) 


1.163(0.193) 


0.615(0.077) 



TABLE II: Dependence of the second order scaling exponent c on the numerical simulation pa- 
rameters: lattice system size L, inverse temperature /?, and nearest neighbour coupling J. The 
exponent c is determined from a fit of the logarithm of the zeros density lnA(i(/)) defined in (|25p to 
the second order scaling function defined in (128p . The value of Ar/ij,, which measures the scale of 
the nonlocal fluctuations of the lattice ensemble, is determined from a fit to (|29|) . 

VII. NONPERTURBATIVE RENORMALIZATION 

Importantly, our new analysis defines the analytic continuation of the lattice partition 
of the continuous-time QMC method as a function of (3 and J. This then implies that the 
second order scaling behaviour that we have identified from the partition functions of our 
analysis corresponds to the zero temperature IR point of the quantum spin chain system 
at 6* = TT, since we have Wick-rotated the spin operator definition via the exact mapping 
j3 — s> ifi. Although we have generated our numerical data at finite temperature, via this new 
analysis, it is therefore possible for us to investigate the zero temperature IR point of the 
quantum spin chain system at 6* = tt directly using finite-temperature lattice ensembles. We 
are able to do this, essentially, by defining the line of constant physics which passes through 
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the point defined by the lattice parameter values in the (3 — J plane and the IR fixed point. 

In {g] and we have recently treated the renormalization of the generic quantum spin 
chain system in the vicinity of this IR fixed point at = tt, using a lattice perturbation 
theory formalism. Recent conformal arguments have suggested that in addition to the usual 
local Sine-Gordon model effective field theory picture of the quantum spin chain in this 
region l]j4|, it is also important to further consider nonlocal quantum effects in this region 
jy generalising these effective field theory treatments to the double Sine-Gordon model 
5]. Our analysis in and suggests that, for quantum fiuctuations of arbitrary size, 
the double Sine-Gordon model should generically require nonperturbative renormalization. 
The analytic continuation properties of the lattice partition we have defined in this article 
essentially define that nonperturbative renormalization program. We quantify the simple 
local fiuctuations of the quantum spin chain via the second order scaling contribution, but 
in addition, we also quantify the relative size of the nonlocal quantum fiuctuations via the 
nonanalytic contribution Atir. Therefore, we can use this scheme to relate any arbitrary 
point in the l3—J coupling plane of the experimental system that this lattice model represents 
to the critical scaling at 6' = tt, and quantify the scale of the quantum fiuctuations in physical 
units. 



VIII. ENTANGLEMENT ENTROPY 



As we discussed in Section 5, the defining relation for the zeroes density of the lattice 
partition function of the quantum spin chain system in f|T9l) is an expression for an entropy. 
Specifically, it is the regularized local limit of the more usual Von Neumann entropy in 
(22), defined via a (^-function renormalization prescription 3J] 35| 36|. To obtain f|T9l) we 
have considered the continuous- J limit of the lattice partition function in (6), in which the 
local transition probabilities Aj become a continuous function of J, given by X{(3,J). We 
have then shown explicitly that this continuous-limit can be reached from a discrete lattice 
system. Importantly, this definition connects with how the Von Neumann entropy in 38|] 37 1 



describes how the subsystem of a quantum spin system is entangled with the remainder of 

these subsystems are defined for our system via the intervals [0, J[ 



the system. From 



and [J, where the full quantum system (universe) is defined via the region [0,f3[. From 
the meromorphic continuation properties rigorously justified in Section VB, it follows that 
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the entanglement entropy of the discrete quantum spin chain lattice system is given by, 



H = -Tr [ A(/5, J) lnA(/3, J) ] (31) 

where the entanglement is defined between the two subsystems of the lattice system 
defined by the bipartite Trotterization of the lattice in (3). 

Several conformal field theory results have been derived for the FSS of the entanglement 



entropy of quantum spin chains 
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44| . The basic result is that the entanglement entropy 
of a quantum spin chain scales like the logarithm of the physical area of the system. This 
area can be parameterised by the physical lattice units of the partition function in ([6]), and 
the second order scaling result we have found is that the scaling of the zeroes density (which 
is directly related to the entanglement entropy via (22)) is directly proportional to the area 
of the lattice in physical units of J, see Table 1. However, the strength of our analysis is 
that we are able to further quantify the nonlocal corrections to this simple scaling pattern 
in the vicinity of the IR fixed point of the quantum spin chain system at ^ = vr, and we are 
able to use then use this scaling result to give precise estimates for the central charge of the 
lattice ensembles we have generated numerically. 

It is known from conformal field theory arguments that the central charge, c, of the general 
quantum spin chain system acquires a very specific value at the point 6 = tt, and that this 
point is related to the unstable renormalization group fiow of the SU(2) Wess-Zumino-Witten 
model and the emergence of a gapless groundstate in half-integer spin chains. According to 
Haldane and Affleck gapped groundstates states are realised by quantum spins chains 
with integer spin values. However, spins chains with half-integer spin values are all gapless, 
and all correspond to the same fixed point of the SU(2) Wess-Zumino-Witten model with the 
same central charge value of c = 1. In the mixed-spin quantum spin model we have treated it 
is known that varying the ratio of the nearest-neighbour spin interaction coupling between 
the spins of different magnitudes (a in (1)) drives the model through the vacuum angle 
6 = TT. This is evidenced by the fact that the finite gap becomes vanishing at specific values 
of this coupling ratio [l^|3][4l| 42|. Varying the ratio couplings in this model, therefore, 
interpolates between these two different physical states of the generalised quantum spin 
chain system: gapped and gapless. Strictly, however, no measurement of the inverse spin 
gap at this gapless point can be truly a divergent function in a finite-temperature lattice 
data set as a function of increasing lattice volume, and this measurement is only strictly 
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a divergent function in these calculations upto the intrinsic IR cutoff scale of the finite 
system. This explicit cutoff scale dependence is quantified in our measurements via Ar//j. 
This nonanalytic term and correction to the area law FSS of the entanglement entropy of 
quantum spin chains has been also identified in Density Matrix Renormalization Group and 
Exact Diagonalization calculations in 45|| 46| 47l |. 

The specific value of the scaling exponent of the entanglement en trop y of quantum spin 
chains is x(c + c)/6 from conformal field theory calculations 38| 4J] 48|, where x is a geo- 
metric factor known as the Euler number which parameterises the topology of the boundary 
between spins, and c is antiholomorphic central charge of the quantum spin chain system. In 
our calculation c = 0, and x = 6. This first value follows because the boundary conditions 
for Euclidean-time are defined as (0,/3], and so only forward propagating modes along the 
Markov chain are defined analytically via (19). This corresponds to the physical problem 



discussed in 
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of the entropy of a conformal system defined by reflections to the right of 
a moving mirror. The second value, for the Euler number, follows from our discussion in 
Section 2. The basic form of the continuous-time method is unaffected by including spin 
vectors of larger magnitude, which can be simply incorporated by a permutation factor |9|] . 
Therefore, even though we have considered a nonuniform quantum spin chain for our analy- 
sis with periodic cells of the form 1 — 1 — 3/2 — 3/2, the geometry of the boundaries between 
two subsystems is unaffected by this choice of mixed-spin model. Nearest-neighbour spin 
sites generically give an Euler number of two 3] and so for the lattice partition function 
of the quantum spin chain defined in (jll) there are essentially three independent nearest- 
neighbour sites, giving a generic x = 6 for this method. Therefore, our conjecture that 
the scaling exponent we have measured describes the second order scaling of the quantum 
spin chain system in the vicinity of the point 6 = n seems to be born out by the numerics. 
The prediction of Affleck and Haldane for the central charge of this IR fixed point in the 
generalised quantum spin chain system is c = 1, and we have measured c = 1 to within 20% 
for all the lattice ensembles we have generated. Discounting the larger volumes at L = 128 
this result is further improved to agreement with unity to within 1% for all data sets. The 
size of the error bars on these measurements are also in good agreement with the results of 
recent DMRG and exact diagonalization studies. 
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IX. SUMMARY 



We have presented a new technique to identify the partition function zeros of quantum 
critical points directly from the information contained in the numerical transfer matrix of 
the continuous-time Quantum Monte Carlo method. We have applied this scheme to a 
mixed spin quantum spin chain model with a nearest neighbour AFM coupling of the form 
1 — 1 — 3/2 — 3/2, and have used the scheme to quantify the renormalization group flow 
of the physical lattice couplings (3 and J of the model in the vicinity of the ^ = vr IR 
fixed point of the quantum spin chain system. Our new technique separates the local and 
nonlocal degrees of freedom of a nonperturbative lattice system analytically, and this has 
allowed us to investigate the renormalizability of the recent proposal that the 6 = ir fixed 
point of quantum spin chains is best described by the double Sine-Gordon model In 
Section 3 we have given a generic treatment of the continuous-time method lattice partition 
function so that our new expansion can readily be generalised to other low-dimensional 
quantum spin systems. In Section 7 we have argued that ( fTOll defines the local limit of the 
entanglement entropy of our quantum spin chain system defined via a bipartite Trotterization 
and, therefore, that we can compare our numerical results with conformal predictions for the 
FSS of the entanglement entropy. We have then done this, and have found good agreement 
between our measured scaling exponents and those predicted by conformal field theory. We 
have identified the correct value of the central charge of c = 1 for the local part of the 
effective field theory description of the quantum spin chain in the vicinity of the 6 = ir fixed 
point of the system predicted by conformal field theory 4J] 38|] , but we have also rigorously 
quantified the nonanalytic correction due to nonlocal fluctuations. Practically this result 
can be used to relate the experimental couplings of a quantum spin chain to the underlying 
conformal symmetries of the system via a program of nonperturbative renormalization. 
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